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1. Introduction 


The study of the cosmic ray fluxes propagation inside meteorites 
is a way to improve our knowledge both on these objects (exposure 
ages, size of irradiated parent bodies, location of samples within 
these bodies..,) and on these fluxes (spectral index, total 
intensity, time variation...). ■ 

In this communication, we present some preliminary results 
obtained from a model built to evaluate the different interactions 
between cosmic rays and meteorites. With this model, we are able to 
compute fluxes as a function of depth inside the meteorite, taking 
into account energy losses by ionization and spallation reactions 
which induce both particle destruction and production. This 
procedure, in which particle fluxes and cross sections are treated 
independently from each other, differs significantly from the thick 
target approach (see Kohman & Bender-1967, Trivedi & Goel-1973) based 
on the measurement of spallation products generated in accelerator 
experiments. 

2. Cosmic rays propagation inside meteorites 

The basic eguation describing the independent of time variation 
with depth of the flux (cm - s~ MeV - ) can be written as : 

d(u)i $j) 

QV % + o n = Qi + — (1) 

' 5E 

where . a is the total destruction cross section of nuclear species i 
. n is the number of target nuclei per unit volume, 

. Q a is the production term of secondary particles, 

. Wj is the energy loss of nucleus i due to excitation and 
ionization of the target medium (u)j= -dE/dx) 

. The streaming factor QV depends on the geometry of the 
meteorite : 

If 9 is the propagation angle, 
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Q V = cos 0 in plane parallel geometry 

5x 

x being the depth inside the meteorite 

n r? a 5 sin 2 0 5 . , , 

Q V = cos 0 — + ~ m spherical geometry, 

dr r 5(cos 0) 

r being the distance from the center. 


3, Asymptotical behaviour of the transport equation at high energy 

Above 1 GeV N“, energy losses and secondary particles 
production become negligible. Therefore, equation (1) reduces to : 


Q V + a n 


0 


The boundary conditions are 
<t> (0, cos 0) = 


(R, cos 0) 


$ 


1.0 


( 2 ) 


(O<0<-) in slab geometry 
2 
7t 

(- _ <0 <it ) in spherical geometry 
2 


(we assume here that the incoming cosmic ray flux $i,o 
isotropic) . 

With such conditions, the solutions of equations (2) are : 


$slab (x,cos9) = exp(- —g.) 

r nr 

“^spherical (r,cos0)= $ 0 expL- a n(r cos 0 + /R ~ 


r sin 


0 ] 


In figure 1, we compare the integrated flux over all directions 
calculated as a function of depth inside a meteorite in the case (i) 

of a plane parallel 

1 1 \ T 1 ' T l 1 t n i i i i i r~a shape and (ii) of a 

spherical shape with 
R= 40 cm. It can be 
seen that the flux is 
always lower if cal- 
culated in the plane 
parallel assumption 
than in spherical 
geometry. The error 
made in treating a 
spherical 40 cm 
radius meteorite as 
an infinite slab 
would only be of a 
factor 1.2 at the 
surface but it would 
(2 £>ecome as high as a 
40 factor 8 in the 
center. 
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Figure 1 

Comparison of depth dependent integrated fluxes 
in spherical geometry' and (2) in slab geometry. Note 
that the surface flux (1) is higher than (2) because 
of outgoing particles. 
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4. Solution of the complete eq uation including secondary particles 

production and energy losses ' ” 

We have adopted the following procedure and assumptions: 

(i) the ionization energy loss dE/dx was calculated from the 
Gloeckler formula (1970), 

(xi) the secondary particles production can be written as : 

Q(E,cos0,r) = £ o(E’) P(E’-*E) ff(E' ,cos0,r)dE' 

where P(E’^E) is the probability to produce a secondary particle of 
energy E from a primary particle of energy E*. 

P 

In this preliminary approach, we used P(E'+E)= — 

E 1 

where P , the average number of outgoing secondaries per spallation 
reaction is an ajustable parameter of our calculation, and the 
secondary spectrum is assumed to be flat. 

The presence of an integral term in the second member of (1) 
requires the use of iterative computations. 

This can be described as follows : 

5o)i <|>(n-1) 

(Q V + a n - — ) 4>( n ) = f° P dE'. 

dE E E' 

Details of the flux computations will be described in Zanda (1985). 

The spallation cross-sections used to deduce nuclei production 
are derived from Silberberg and Tsao (1972a, 1972b, 1977, 1979). 

The energetic particles may have two origins : 

(i) the galactic cosmic rays which are modulated by the solar 
activity (we derived our modulated spectrum out of Proteroe, Ormes 
and Comstock -1981), 

(ii) the solar cosmic ays. Since this solar source consists in 
particles with E<100 MeV, it only affects the first few centimeters 
of the meteorite and will not be considered here. 

5 . Preliminary results 

Because of the existing experimental data, we chose to first 
test our model on the iron meteorite Grant. Computations were made 
for a meteorite having the same chemical composition and various 
radius (40 cm, 35 cm and 30 cm) admiting a secondary particle yield 
per interaction p=0.75. Figure 2 shows a comparison of our results 
with measurements by Signer and Nier (1960) for Ne21 . It can be seen 
that there is a good agreement between calculation results and 
experimental datas for a radius of 35 cm which is a little less than 
the pre-atmospheric radius that these authors derived of their 
calculation. 
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Figure 2 

Comparison of computations results for various radius (single curves) 
and experimental datas from Signer and Nier (dotted curve). 


6 . Conclusion 

The results obtained with the model presented here are sufficiently 
encouraging to pursue and improve that type of studies, in particular 
by finding a proper treatment for the outstanding problem of 
secondary particles production. Unfortunately, taking into account 
the actual geometry of a meteoritical sample at the time it suffers 
the bulk of the cosmic ray irradiation will still remind a 
difficulty. 

We wish to thank G. Malinie for helpful discussions throughout the 
preparation of this paper. 
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